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ABSTRACT 

We examine how the statistics of the quadrupoles of (projected) cluster masses can discriminate 
between fiat cold dark matter (CDM) universes with or without a cosmological constant term. Even in 
the era of high precision cosmology that cosmic microwave background experiments should open soon, it is 
important to devise self consistency tests of cosmogonic theories tuned at the matter radiation decoupling 
epoch using data from the non-linear evolved universe. We build cluster catalogs from two large volume 
simulations of a "tilted" CDM model and a ACDM model with cosmic density parameter Q, m = 0.35 and 
cosmological constant contribution Ha = 0.65. From the projected mass distribution of the clusters we 
work out the quadrupoles Q and examine their dependence on cluster mass and the cosmological model. 
We find that TCDM clusters have systematically larger quadrupoles than their ACDM counterpart. The 
effect is mass dependent: massive clusters (M> 10 15 /i _1 Af Q ) have quadrupoles differing by more than 

30% in the two models, while for M^S 4 x 1O 14 /i _1 M0 the difference rapidly drops to ~ 1%. Performing 
a K-S test of the Q distributions, we estimate that using just the 15 most massive clusters in the 
simulation volume (360 /i _1 Mpc a side) we can discriminate between TCDM and ACDM at a confidence 
level better than 99.9%. In the volume probed by exhisting observations, there are potentially several 
hundred clusters with masses above the threshold for which the differences in the quadrupoles become 
relevant. Should weak lensing data become available for this whole set, a quadrupole analysis may be 
expected to discriminate among different values of A. 

PACS: 95.35; 98.80; 98.65.Cw 

Subject headings: cosmology: theory - dark matter - large-scale structure of the Universe - galaxies: 
clusters - galaxies: halos - methods: numerical 



1. INTRODUCTION 



Clusters of galaxies have been systematically studied in 
recent years, using both optical and X-ray data. In princi- 
ple, it is possible to obtain from their properties stringent 



function (e.g. 


Eke, Cole k Frenk 1996, Eke, Cole, Front 


& Henry 1998, 


Cole ct. al 1997, Viana & Liddlc 1996|, Mo. 


Jing & White 199q, 


Borgani ct al. 1997|, Liirardi et al. 


1998, Postman 1998], 


Gardini et al. 1999|). 



ters have complex morphologies with strong evidence of 
substructure (e.g. Geller & Beers 1982; Dressier & Schect- 
man 1988; West & Bothun 1990; Forman & Jones 1990; 
Bird 1994; West, Jones & Forman 1995; Bardelli et 
al.1998; Solanes et al.1999). In parallel to the develop- 
ment of these observations of the visible (barionic) matter 
in clusters, the study of the distortions of the images of 
distant galaxies gravitationally lensed by cluster potentials 
have made it possible to systematically map the distribu- 
tion of the mass, including the dynamically dominant dark 
matter component (e.g. Tyson, Valdes & Wrenk 1990; 
Kaiser & Squires 1993; Seitz & Schneider 1995; for ongo- 
ing surveys, e.g. Clowe et al.2000, Dahle 2000, Graham 
et al.2000, Wittman et al.2000; for recent reviews Hattori, 
Kneib & Makino 1999, and Kaiser 1999). 

In hierarchical bottom up cosmological models, like the 
cold dark matter model (CDM) and its variants, galaxy 



clusters are the largest bound structures to form (White 
& Rees 1978; Davis et al.1985). Their properties can 
retain the signature of the cosmological parameters more 
easily than smaller older structures. 

Indeed, in different cosmologies clusters can have sig- 
nificantly different formation histories and this may affect 
their typical morphologies. This fact was first examinated 
using analytical methods by Richstone, Loeb & Turner 
(1992), Bartelmann, Ehlers & Schneider (1993) and Lacey 
& Cole (1993), and, using numerical methods, by Evrard 
et al.(1994) Mohr et al.(1995) and Wilson, Cole & Frenk 
(1996b; hereafter WCF96). The latter authors studied the 
dependence on cosmology of the quadrupole of a cluster's 
projected mass by simulating weak gravitational lensing 
in artificial clusters (Wilson, Cole & Frenk 1996a) grown 
in numerical simulations of the CDM model with differ- 
ent values of today's matter density parameter fi m and 
cosmological constant A. The relation between the pro- 
jected mass distribution of clusters and the cosmological 
model has also been recently explored in the context of 
hydrodynamical simulations by Valdarnini, Ghizzardi & 
Bonometto (1999). 

In a low-density universe without cosmological constant, 
the linear growth of density fluctuations ceases after a red- 
shift z ~ (l/fi m ) - 1 (e.g. Peebles 1980, sees 11 and 13). 
If the model is spatially flat thanks to a non-vanishing 
cosmological constant, the growth of fluctuations freezes 
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when the vacuum energy density begins to dominate the 
matter density p m , that is, as p m oc (1 + z) 3 , at a redshift 
z ~ (n^/Qm) 1 / 3 — 1. Hence, in low-density universes, clus- 
ters on average form at moderately high redshifts (z> 0.5) 
and subsequently accrete little material; on the contrary, 
in the standard fi m = 1 CDM universe, structure for- 
mation occurs continuously with rich galaxy clusters also 
forming in very recent epochs (z%> 0.3). WFC96's analysis 
showed that the statistics of the quadrupoles of the pro- 
jected mass distribution of clusters is sensitive to the value 
of O m , quite independently of the value of A. However, 
other than Q m = 1, they considered only a rather extreme 
value fi m = 0.2. Also, they examined only a small sample 
of (eight) simulated clusters rather than building a catalog 
from a large volume simulation and could not examine the 
importance of the effect as a function of cluster mass. 

In this paper we extend the study of how the statis- 
tics of cluster mass quadrupoles is sensitive to the cos- 
mological parameters building cluster catalogs from large 
scale A^-body simulations. In agreement with recent 
data on the angular power spectrum of CMB anisotropy 
(BOOMERANG-98, de Bernardis et al.2000, Lange et 
al.2000; MAXIMA- 1, Hanany et al.2000, Balbi et al.2000; 
see also White, Scott & Pierpaoli 2000, Jaffe et al.2000), 
we consider two flat universes (fi m + ft a = 1). One is a 
standard cold dark matter model with a "tilt" in the pri- 
mordial spectral i ndex (n = 0.8); t his simulation is the 



same analyzed by Gardini et al 
has i^TD 



1999 



The other model 
0.35 and cosmological constant Ha = 0.65 



(ACDM), which fits CMB quadrupole and cluster abun- 
dance data and is compatible with the data on high- 
redshift Type la supe rnovae which suggest an accelerat- 
ing cosmic expansion ( Perlmutter et al. 1998] , Riess et al 
1998Q). The TCDM model is ruled out by combined CMB 
and LSS data (e.g. ), but roughly reproduces the CMB 
quadrupole a nd the observed abu ndance of rich galaxy 



clusters (e.g. |Gardini et al. 1999 ). This is sufficient for 



the purposes of this analysis, since we want to test ACDM 
against an example of a standard (Q = l) CDM cosmology. 

In the next few years, starting from the results of MAP 
(e.g. Page 2000), and also thanks to the Planck experi- 
ment (e.g. De Zotti et al.2000), a detailed reconstruction 
of the angular spectru m of CMB will be performed. Sev- 



eral a uthors (see, e .g., Jungman et al.1996 , Kamionkowski 
& Koi iowsky 1999 ) have shown that Q m and A, as well as 



the spectral index n of primeval density fluctuation, will 
then be safely determined with unprecedented precision. 
However, not all degeneracy in the parameter space will be 
removed, unless suitable a priori conditions are set. For 
instance, small components of hot dark matter, consistent 
with current limits on neutrino masses coming from neu- 
trino mixing experiments, could introduce residual uncer- 
tainties. A problem could also arise if the vacuum energy, 
accounted for by the cosmological constant A, is not re- 
ally constant in time. Indeed, interpreting the data on 
high-redshift supernovae in terms of a cosmological con- 
stant may be simplistic. Such data seem to require that 
the energy density is dominated by a negative-pressure 
component (dark energy); but a constant vacuum energy 
is just one possibility. A popular alternative scenario is 
quintessence (Caldwell et al.1998), which under certain 
conditions can correspond to a time-dependent cosmolog- 



ical "constant" (e.g. Crooks et al.2000). 

This case is potentially dangerous, as the very equation 
of state of dark energy could vary with time. If this occurs 
and the variation is significant, a large part of the informa- 
tion contained in the CMB angular spectrum could have 
to be "spent" to reconstruct the time dependence of the 
components of the stress-energy tensor of dark energy. Al- 
though, in this paper, we do not consider models with time 
dependent A, the results obtained with the technique pre- 
sented here could become particularly significant for such 
a case, providing information complementary to CMB on 
the time dependence of A. 

However, independently of the occurrence of such a del- 
icate case, it should be remarked that the imprints left 
by density fluctuations on the CMB are memories of a dis- 
tant past of the Universe. Data from the local universe can 
be used to provide self-consistency tests for structure for- 
mation theories tuned at the matter-radiation decoupling 
epoch which is probed by CMB observations. 

The plan of the paper is as follows. In the next sec- 
tion we describe the simulations used. In section 3, we 
illustrate how we identify the clusters and measure their 
properties. In section 4 we work out the statistics of clus- 
ter mass quadrupoles and present our results. Finally, we 
summarize and discuss our conclusions in section 5. 

2. THE SIMULATIONS 

We carried out two simulations using the parallel N- 
body code of Gardini et al.(1999), which was developed 
from the serial public AP3M code of Couchman (1991) ex- 
tended to different cosmological models and different mass 
particle sets. The first simulation was a "tilted" Einstein- 
de Sitter model (hereafter TCDM), while the second sim- 
ulation was a ACDM model, i.e. a flat CDM universe 
with non-zero cosmological constant. The parameters of 
the model are reported in the Table. Both these models 
roughly yield the co rrect abundance of rich clusters (see 
Gardini et al. 1999[ ). The TCDM simulat ion is the same 
already considered by Gardini et al. 1999; the normaliza- 
tion of the run has been rescaled to yield as = 0.55 at the 
final epoch, instead of cj 8 = 0.61, so that the abundances 
of rich clusters are the same for both models at the final 
epoch. 

The simulated volumes are 360/i _1 Mpc cubic boxes. 
CDM+baryons are represented by 180 3 particles, whose 
individual mass is 2.22 • lO 12 ^" 1 M for TCDM and 0.777 • 
lO^h^M® for ACDM . We use a 256 3 grid to compute 
the FFTs needed to evaluate the long range contribution 
to the force (PM) and we allow for mesh refinement where 
the particle density attains or exceeds ~ 30 times the mean 
value. The starting redshifts are Zi n = 10 for TCDM and 
z = 20 for ACDM . The particle sampling of the den- 
sity field is obtained applying the ZcPdovich ap proxima- 
tion flZcl'dovich 1970| , |Doroshkcvich et al. 1980| ) starting 
from a regular grid. We adopt the same random phases in 
both simulated models. 

The comoving force resolution is given by the softening 
length, rj ~ 112/i _1 kpc. The force F is evaluated con- 
sidering each particle as a smoothed distribution of mass, 
with shape p(r) = (48/7T77 4 )(t7/2 — r) for r < rj/2 (thi s 



is the so-called S2 shape, Hockney & Eastwood 1981) 



It behaves as F oc 1/r 2 when r > rj. Since the soften- 
ing of the force is usually referred to a Plummer shape, 
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F oc r/(r 2 + e 2 ) 3 / 2 , we used a least x 2 test to estabilish 
the best approximation between the forces generated by 
the two different shapes. The minimum % 2 occurs when 
?7 = 2.768e. In our case, this corresponds to a Plummer 
equivalent softening e p i ~ 40.6h~~ kpc; we will use this 
latter value of the softening as our nominal force resolu- 
tion. Our comoving force and mass resolutions approach 
the limits of the computational resources of the machine 
we used (the HP Exemplar SPP2000 X Class processor of 
the CILEA consortium at Segrate-Milan) . 

The number of steps were 1000 equal p-time steps (the 
time parameter is p oc a 2 / 3 , where a is the expansion fac- 
tor). Such step choices were dictated by two criteria: (i) 
Energy conserv ation. According to L ayzer Irvine equa- 
tions (see, e.g., Efstathiou ct al. 1985 ), it had an overall 
violation < 3 % for TCDM and ACDM . (ii) Cole et al. 
(1997) requirements, that the rms displacement of parti- 
cles in a step is less than ry/4 and the fastest particle has 
a displacement smaller than rj were never violated. 



Nf corresponding to a mass threshold 3.0 x 10 13 /i _1 M s . 
Above this mass threshold there are ~ 10000 clusters in 
the samples. For the analyses, we only use clusters more 
massive than 4.2 x 10 14 /i _1 M© (see § 5); there are ~ 300 
of them. 

4. CLUSTER QUADRUPOLES 

For each cluster identified as above, we compute the 
quadrupole of the cluster particle distribution projected 
along a randomly oriented direction through the simula- 
tion box. 

If Xi and yi are the coordinates of the i-th cluster par- 
ticle when projected onto the plane perpendicular to the 
axis chosen (with the origin at the cluster's center), the 
quadrupole of the projected particle distribution is defined 

as: 



Q 



(g/AW) 1/2 



(1) 



3. CLUSTER IDENTIFICATION 

One of the main aims of the simulations is that of ob- 
taining a large set of model clusters for each cosmological 
model at different redshifts. This will enable us to study 
cluster evolution and to create mock cluster catalogs. Here 
we shall report some basic results and general properties 
of the clusters selected in the simulation outputs. 

The clusters we consider were found using a spherical 
overdensity (SO) algorithm, yielding the cluster locations, 
the radii R s inside which a density contrast S cr — S V i r is 
attained and the total mass M of the particles within R s . 
As "virial" overdensities A vir for TCDM and ACDM we 
choose the "standard" values 178 and 110 suggested by the 
spherical infall mode l; we use the relation A V j r — 17 8f^ 45 
for fi m +fi A = 1 (e.g. |Eke, Cole, Frenk fc Henry 1998Q . The 
SO procedure has the benefit of providing clusters which 
appear to satisfy a sensible virialization requirement (a 
given overdensity in a sphere). 

Let us now describe our SO procedure implementation. 
As a first step, candidate clusters are located using a stan- 
dard FoF algorithm, with linking length A = cf>xd (here d is 
the average particle-particle separation), yielding groups 
with more than Nf particles. We then perform the follow- 
ing operations: (i) we find the center-of-mass Cm of each 
group and (ii) we determine the radius R g , inside which the 
density contrast is 5 cr (all particles are included, not only 
those initially found by FoF). In general, the new center- 
of-mass is not Cm- The operations (i) and (ii), define a 
new particle group, on which the same operations (i) and 
(ii) can be repeated. The procedure is iterated until we 
converge onto a stable particle set. If, at some stage, the 
group contains less than Nf particles, we discard it. The 
final R g is R s . It may happen that a particle is a potential 
member of two groups; in this case the procedure assigns 
it to the more massive one. This has the consequence that, 
sometimes, more massive groups swallow smaller ones, 
producing a slight decrease of the total number of clusters, 
over all mass scales (see Gardini et al.1999). Gardini et 
al.1999 describe the SO algorithm and its outputs in detail 
and compare it with those found using other group identi- 
fication algorithms and previous work (e.g. Govcrnato ct) 
al . 1999] ). In this work SO was started setting <p = 0.2 and 



where 



4(^2/^ 



(2) 



The sum extends over the N V i r particles contained within 

Rvir • 

For the analyses presented in the next section, we have 
repeated the calculations using three random directions to 
compute the quadrupoles and verified that our conclusions 
do not depend on the projection chosen. The results are 
presented for one projection axis. 

The scales resolved by our simulations are significantly 
larger than those resolved in recent high resolution works 
(Brainerd et al. 1997, Moore et al.1998, Ghigna etal.1998, 
Klypin et al.1999, Okamoto & Habe 1999). Therefore, 
it would be hard to trace in our clusters the rich mass 
substructure found by these authors. It is however un- 
likely that this affects the evaluation of a global property 
such as the quadrupole. A rough estimate of the resolu- 
tion needed to obtain accurate quadrupole measures can 
be performed by comparing the linear scale resolved (r res ) 
with ttR/2 (R is the distance from the center of the clus- 
ter considered) . Only the contribution to the quadrupoles 
coming from substructure in a central region of radius 
not exceeding a few hundred kpe's should therefore be af- 
fected by lack of resolution. Let us consider a cluster of 
radius R and the contribution to its quadrupole coming 
from a central region of projected radius r. The ratio be- 
tween the volumes of the central region and the whole 
cluster is ~ (3/2)(r/i?) 2 . With R ~ 1.5/i _1 Mpc and 
r ~ 150 ft. _1 kpc, such ratio is ~ 1-2 %. The mass profile is 
not flat and this increases the weight of the central region; 
however, this bias should not be important because the 
structure of the central regions should be mostly shaped 
by the internal dynamics of the cluster, with little memory 
of late tidal actions and/or dishomogeneities due to sec- 
ondary infall, which most likely cause the differences in the 
global anisotropy of the clusters in the different models. 

Besides possible numerical effects, there is another lim- 
itation in our analysis. The noise that would affect the 
estimates of Q from the mass distribution of real clusters 
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has not been taken into account (e.g. from projection ef- 
fects and from the limitations of the techniques used to 
reconstruct the mass from the observed gravitationally in- 
duced shear; see Wilson, Cole & Frenk 1996a). We plan 
to improve the analysis along these lines in future work. 
However note that the analysis of Rcblinsky & Bartelmann 
(1999), using large scale cosmological simulations, indi- 
cates that projection effects in weak-lensing-mass selected 
cluster samples should be small. 



5. STATISTICAL ANALYSES 

For both cosmological models, we have examined the 
distribution of the quadrupoles. As pointed out in the 
Introduction, these estimators of cluster morphologies are 
expected to probe the different formation histories of clus- 
ters in different cosmologies. In both cases we select the 
clusters more massive than 4.2 • 10 14 M Q ; the numbers of 
clusters above threshold are the same for the two models. 

The histogram in Figure ^ shows the overall distribu- 
tion of the quadrupoles measured. TCDM clusters have 
clearly larger values of Q than their ACDM counterparts, 
in agreement with the expectation of ongoing infall and ac- 
cretion of present clusters in the TCDM model. This effect 
is more prominent for very massive clusters, as it is shown 
in Figure | The fi gure shows the trend of Qvs.M V { r . 
Each point is obtained ranging the clusters in order of in- 
creasing mass and averaging over 100 neighbours in the 
lists; the errorbars correspond to one standard deviation 
in the distribution of quadrupoles. 
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Fig. 1. — Histogram of the distribution of the quadrupoles for the 
whole set of clusters more massive than 4.2 X 1O 14 M (On the ver- 
tical axis we simply plot the number of clusters above threshold in 
each bin). TCDM clusters have systematically larger quadrupoles 
than their ACDM counterparts. 
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Fig. 2. — Average Q vs. virial mass M v i r . The value of Q 
is obtained ranging the clusters in order of increasing masses and 
averaging over 100 neighbours (i.e. they are binned into bins con- 
taining 100 clusters each): the value M v i r reported is the average 
of the 100 values. 
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Fig. 3. — Scatter plot of Q vs. M„; r for rich clusters. 

The average quadrupoles of TCDM and ACDM clus- 
ters are sistematically different at all masses; the most 
interesting result however is the strong difference found 
at the high mass end. Note that this effect could not be 
detected by WCF96, because of the small number of simu- 
lated clusters in their analysis and, most importantly, be- 
cause they built a sample of essentially equal mass clusters 
(they rescaled the particle masses of their runs so as to ob- 
tain same mass objects). Figure |^ is a scatter plot of Q vs 
M V i r for clusters more massive than 6.3- 1O 14 M0, the mass 
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threshold for which the average quadrupoles of TCDM and 
ACDM clusters start differing by several percentual points 
(this threshold yields f05 clusters for ACDM and 81 for 
TCDM). 

How significant is the difference between the two dis- 
tributions ? Could it be detected through observations 
of real clusters ? Observations of clusters have made im- 
pressive progress in the past decade. As we have reviewed 
in the Introduction, maps of the mass distribution itself 
reconstructed from observations of gravitational lensing 
effects can be systematically built. On the observational 
side, much progress is a consequence of the development of 
large-format CCD cameras which can cover fields of ~ 0.5 
degrees on a side on 4-meter class telescopes. On the 
theoretical side, different inversion algorithms have been 
developed to systematically exploit the data on weak grav- 
itational shear, accounting for the effects of realistic PSFs 
and optimally weighing background galaxies (see the re- 
cent reviews by Mellier, 1999, and Bartelmann & Schnei- 
der, 2000). Compilations of existing data would allow sam- 
ples of ~ 40 clusters to be examined (e.g. Mellier 2000). 
However, such samples are hardly homogeneous and the 
analysis of them would require a careful examination of 
systematic effects. Only very recently an effort to per- 
form a weak lensing study of a large well defined sample 
has been started using the NOT and the University of 
Hawaii 2.24-m telescope (Dahle 2000). On the contrary, 
inspecting a few very massive clusters is a relatively easy 
task. This is particularly welcome for the purposes of the 
present analysis, given that the most prominent differences 
between the models appear at the high mass end. Let us 
focus on the behaviour of, say, the first 15 most massive 
clusters in both samples. As seen in Figure | TCDM clus- 
ters are clearly shifted towards high values of Q. 




Clusters Number 

Fig. 4. — Kolmogorov-Smirnov probability that the Q measured 
for the N most massive clusters in the simulation volumes are drawn 
from the same theoretical distribution. The (red) dots without 
errorbars refer to clusters selected in the whole simulation boxes 
(360/i _1 Mpc a side). The (blue) dots with errorbars are averages 
over the results obtained for (eight) half-size subboxes (180 h — Mpc 
a side); the errobars correspond to one standard deviation. 



In order to quantify the significance of this difference, we 
perform a K-S test on the distribution of the quadrupoles. 
Figure |^ gives the probability that the Q measured for the 
N most massive clusters in ACDM and TCDM are drawn 
from the same underlying distribution (red dots, without 
errorbars). The cluster number N is varied between 4 and 
15. Even with only 10 clusters the probability that the 
difference between the Q measured for the ACDM clusters 
and their TCDM counterparts is just a statistical fluke is 
below 1%; using 15 clusters, it is well below 1/1000. Ex- 
tending the analysis to all the clusters more massive than 
6.3 • 10 14 M Q (about twice the mass of the Virgo cluster) 
yields a KS probability of a mere 10~ 32 . 

The simulated clusters live in a cubic volume of 
360 /i -1 Mpc a side. According to the result above, picking 
~ 10 massive clusters (M V i r ~ 10 15 M Q ) in such a volume 
is sufficient to discriminate between the two models at the 
1% confidence limit. Of course, the result will depend on 
the volume explored, both because 10 15 M Q clusters are 
rare objects and because of cosmic variance. Ideally, differ- 
ent simulations would be welcome to thouroughly examine 
this issue. As a tentative test, we have divided the simu- 
lation boxes in eight subboxes of 180 A. Mpc a side; for 
each pair of corresponding ACDM and TCDM subboxes, 
we have then repeated the same analyses as for the parent 
boxes, i.e. we have selected the 15 most massive clusters 
and estimated the KS probabilities that their Q distribu- 
tions are drawn from the same underlying model. The 
results are shown in figure ^ as blue dots with errorbars 
(l-cr). The average KS probability increases markedly; 
however the result is dominated by the outcome of the KS 
analysis in one of the subboxes, in which the most mas- 
sive clusters have very similar distributions, whereas in 
the other subboxes the KS test yields probabilities only 
slighlty higher than for the whole box. This indicates that 
at least a volume > 200/i _1 Mpc is necessary to obtain 
robust quadrupole statistics able to discriminate between 
the two cosmologies just using a few tens of massive clus- 
ters. A set of large volume simulations would be welcome 
to examine this issue further. 

6. DISCUSSION AND CONCLUSIONS 

Let us first summarize the analyses reported in this pa- 
per. We have used two large volume simulations of a 
ACDM and a TCDM universe, which were started from 
the same initial random numbers so that the structure 
forming in each of them are directly comparable. The 
simulations span a cubic volume of 360 Mpc a side. The 
numerical parameters of the simulations and the linear fea- 
tures of the cosmological models are examined in § 2. 

In each simulated volume we have identified the mas- 
sive virialized systems that would correspond (in mass) 
to observed galaxy clusters. We have then measured the 
quadrupole Q of the projected mass distribution of each 
simulated cluster (we extend our analysis to objects more 
massive than a threshold = 4 x 10 14 M Q ). 

Studying the statistics of the quadrupoles, we have ob- 
tained the following results: 

• The distributions of the quadrupoles of TCDM 
and ACDM clusters differ significantly, with TCDM 
clusters having systematically larger quadrupoles in 
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agreement with the linear theory expectation that in 
a high density universe clusters keep accreting ma- 
terial till the present epoch while in a low density 
universe late accretion and mergers are suppressed. 
This confirms previous numerical results by Evrard 
et al.(1993) Mohr et al.(1995) and Wilson, Cole & 
Frenk (1996b). 

• At variance with the previous authors which focussed 
on a small set of (at that time) high resolution sim- 
ulations, we have built catalogs of simulated clus- 
ters. This allows us to study the dependence of the 
quadrupoles on cluster mass. We find that the trends 
of Q vs. M V i r of the two models are markedly dif- 
ferent at the high mass end. For M V i r > 8 • 10 14 M Q 
(about twice the mass of the Virgo cluster), the av- 
erage quadrupoles of ACDM clusters drop to values 
30% lower than those of their TCDM counterparts 
of similar mass. 

• We have examined how many massive clusters 
are needed to significatly distinguish between the 
quadrupole distributions of TCDM and ACDM clus- 
ters performing a KS test. We find that, using 
just the 10 most massive clusters of the samples 
(> 1O 15 M in the simulation volume of 360 ft, _1 Mpc 
a side), the probability that the Q measured are 
drawn from the same model distribution is less than 
1%, and decreases to well below 1/1000 if the 15 
most massive clusters are employed. 

• The significance of the test depends on the volume 
span by the cluster sample and on the number of 
most-massive clusters which can be taken into ac- 
count. If such number is restricted to O(10), we 
estimate that a volume > 200/i _1 Mpc is necessary 
to obtain robust quadrupole statistics able to dis- 
criminate between the two cosmologies. Of course, 
increasing the number of clusters (that is lowering 
the mass/richness threshold for a cluster to be con- 
sidered) implies more and more careful observational 
work. It should also be noted that, if the mass 
threshold is pushed down to 4 x 10 14 /i _1 M©, no 
substantial advantage would result from adding more 



clusters to the sample, as, at such masses, the sig- 
nal tends to quickly fade. For example, in a box 
with side ~ 180 h~ 1 Mpc, we find ~ 25 clusters above 
such threshold. In a volume ~ 10 7 -10 8 /i~ 3 Mpc 3 , of 
the order of what is observationally inspected, the 
number of significant clusters approaches a total of 
~ 200. A set of new simulations could be used to de- 
termine how effectively different values of A can be 
discriminated, should this whole sample of potential 
data be exploited. From our analysis, we speculate 
that it could allow us to go beyond discriminating 
between A = and A / 0. 

These results confirm that the statistics of cluster 
quadrupoles is sensitive to the cosmological model. In par- 
ticular this statistics could do more than distinguishing be- 
tween the currently popular ACDM model and a standard 
CDM cosmology with zero cosmological constant. Even 
in the era of high precision cosmology that CMB observa- 
tions will be hopefully opening soon, the statistical anal- 
ysis presented here might be quite significant. First of 
all, it conveys information coming from our own or nearby 
epochs, allowing a self consistency test of theories tuned 
at the matter-radiation decoupling epoch. Furthermore, 
in the case of complex models of dark energy as discussed 
in the Introduction, it can complement CMB observations, 
helping to remove residual degenerations in the determi- 
nation of the cosmological parameters. Of course, us- 
ing quadrupole statistics requires, from the observational 
point of view, a strong systematic effort to expand the 
available data set on weak lensing from galaxy clusters, 
and, from the theoretical side, a careful examination of 
the biases affecting the quadrupole estimates of real clus- 
ters and their selection (e.g. due to projection effects or 
the limitations of cluster identification and mass recon- 
struction algorithms; see e.g. Wilson, Cole & Frenk 1996a 
and Reblinsky & Bartelmann 1999). 

We thank the consortium CILEA for allowing us to run 
a former trial simulation free of charge and a particular 
thank is due to Giampaolo Bottoni of CILEA, for his ex- 
pert technical assistance. 
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Constraining A using Cluster Quadrupoles 





TCDM 


ACDM 




1 


0.65 


A 





0.35 


n b • io 2 


6 


6 


n 


0.8 


1.05 


h 


0.5 


0.65 


QpS,rms/ 


15.6 


20.05 




0.55 


1.08 


r 


0.32 


0.19 


Ar c/ (PS; <5 C = 1.69) 


4.2 


4.2 



Table 1 

Parameters of the mod els. All parameter s listed are either input parameters or quantities worked out from 
the linear theory (see gardini et al. 1999). the normalization to cobe quadrupole was deliberately kept at 
the ~ 3(t lower limit, in order to leave some room to the contribution of tensor modes, while being consistent 

with the data. the expected interval for n ct is 4-6. 



